METHOD AND APPARATUS FOR CALIBRATING AND MEASURING 
ARTERIAL COMPLIANCE AND STROKE VOLUME 



5 Related Applications 

This is a divisional of patent application Ser. No 10/21 1,857 filed Aug. 1, 2002. 

Field of the Invention 

This invention relates to the field of medical diagnosis, and more specifically, to a 
10 method and apparatus for noninvasive measurement of arterial compliance. 

Background of the Invention 

High blood pressure and hardening of the arteries can lead to coronary problems. 
Obtaining various measurements of the vascular system, including compliance of large 

15 and small vessels, and systemic resistance, provides physicians with information useful in 
diagnosing and treating early stages of coronary disease. 

U.S. Patent No. 6,017,313 ("the 313 patent") to Christopher W. Bratteli et al. 
(incorporated herein by reference) discloses an apparatus and method for blood pressure 
pulse waveform contour analysis. Methods of processing an arterial blood pressure 

20 waveform to extract clinically useful information on the state of the cardiovascular 
system are disclosed. In order to obtain the parameters of the modified Windkessel 
model, the diastolic portion of a subject's blood pressure waveform is scanned over a 
plurality of ranges and the range that produces the best fit of data and lowest error 
estimates are selected. 

25 U.S. Patent No. 5,21 1,177 ("the '177 patent") to Charles F. Chesney et al. 

(incorporated herein by reference) discloses method and apparatus for measuring 
properties of the human vasculature using an electrical analog model of vascular 
impedance. These properties include the compliance of large and small vessels, and 
systemic resistance. These measurements and others obtained from the model can, in 

30 turn, be used to diagnose states of health or disease, and to assess the effectiveness of 
treatment regimes. For example, see Finkelstein, S.M., Collins, V.R., Cohn, J.N., 
Arterial vascular compliance response to vasodilators by Fourier and pulse contour 
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analysis, Hypertension 1988:12:380-387, the entire disclosure of which is incorporated 
herein by reference. 

U.S. Patent No. 6,159,166 ("the '166 patent") to Charles F. Chesney et al., and 
entitled SENSOR AND METHOD FOR SENSING ARTERIAL PULSE PRESSURE 
5 (incorporated herein by reference) discloses a method and apparatus useful for tonometric 
measuring of an arterial pulse pressure waveform. U.S. Patent No. 6,331,161 to Charles 
F. Chesney et al., and entitled METHOD AND APPARATUS FOR FABRICATING A 
PRESSURE-WAVE SENSOR WITH A LEVELING SUPPORT ELEMENT 
(incorporated herein by reference) discloses a method and apparatus of a sensor useful for 

10 tonometric measuring of an arterial pulse pressure waveform. U.S. Patent No. 6,132,383 
to Charles F. Chesney et al., and entitled APPARATUS AND METHOD FOR 
HOLDING AND POSITIONING AN ARTERIAL PULSE PRESSURE SENSOR 
(incorporated herein by reference) discloses a method and a wrist brace and sensor holder 
useful with any of the above sensors for tonometric measuring of an arterial pulse 

1 5 pressure waveform. 

U.S. Patent No. 5,241,966 ("the '966 patent") to Stanley M. Finkelstein; et al. 
issued Sept. 7, 1993 (incorporated herein by reference) discloses an apparatus for 
measuring stroke volume/cardiac output that includes a transducer for measuring arterial 
blood pressure waveform, a digitizer for digitizing the analog signal generated by the 

20 transducer and a digital signal processor for determining ejection time and heart rate. 

Processor circuitry determines cardiac output using the ejection time, heart rate, the body 
surface area and age of the patient. 

Compliance is a fundamental property of any pressure- volume system such as the 
arteries. In the arteries, reduced compliance has been offered as a mechanism whereby 

25 the work load on the heart is increased and myocardial perfusion is decreased, thus 

leading to cardiac disease. It is also thought that disease of the arteries themselves may 
be evaluated by measurement of arterial compliance. Therefore, arterial compliance is an 
important cardiovascular parameter. 

What is needed is improved noninvasive measurement of arterial compliance. 
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Summary of the Invention 

The present invention provides improved noninvasive measurement of arterial 
compliance using a combination of noninvasive arterial tonometry and noninvasive cuff 
oscillometry. Described are a number of improved approaches to estimating systemic 
5 vascular resistance and/or compliance. In one embodiment, a computer-controlled 
pneumatic cuff having a sensor coupled to sense pressure and pressure variations in the 
cuff is used to obtain oscillometric signals (e.g., from an upper arm of a patient), and 
simultaneously a contact-pressure sensor (e.g., some type of microphone) is used to 
obtain tonometric signals (e.g., from the contralateral wrist of the patient and, e.g., a 
10 transducer placed on the radial artery), and the oscillometric signals and tonometric 
signals, which are correlated to one another in time, are combined to calibrate one 
another and/or obtain information about the patient not available conventionally. In some 
embodiments, the oscillometric measurement is used to calibrate the tonometric 
measurement. 

1 5 One aspect of the present invention provides a method to compute arterial 

compliance from oscillometric data (2 nd embodiment listed below). This provides a 
compliance measurement that is distinct from the compliance measurement described in 
the '177 patent. 

Another aspect of the present invention provides a method for improving pressure 
20 calibration from oscillometric data (1 st and 3 rd embodiments listed below). This can be 
used to improve the above measurement, to improve the ! 177 patent's described 
compliance measurement, or to improve almost any measurement that uses arterial 
pressures. 

Yet another aspect of the present invention provides a set of methods for building 
25 models that compute flow from pressure data (with or without oscillometric data). This 
can be used in the same way as the cardiac output model described in the 313 patent (i.e., 
it can be used to improve the compliance estimate). 

The present invention provides a number of improvements or refinements to the 
approaches to waveform analysis set forth in the '177 and the f 313 patents. These 
30 improvements include an accurate calibration method and apparatus for using an 

oscillometric signal to calibrate the pressures of tonometric signals in the contralateral 
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arterial site. These improvements also include a simple uncorrected volume arterial 
compliance curve that is obtained by plotting relative arterial volume under the cuff 
against brachial arterial transmural pressure. Further improvements account for a bias 
(i.e., an overestimation of the in vivo compliance curve for the arterial segment due to 
5 shear stresses at the ends of the cuff) by providing a correction for the transmural 
pressure that is based on stress-strain properties of the upper arm. 

In some embodiments, the method of the '313 patent is improved with use of the 
SVR model of the present invention. 

10 Brief Description of the Drawings 

FIG. 1 A illustrates an example of a vascular-compliance-determining device 100 

according to the present invention. 
FIG. IB illustrates schematically details of one embodiment of oscillometric sensor 150. 
FIG. 1C illustrates schematically details of one embodiment of cuff 158 over an artery. 
15 FIG. 2 A illustrates an oscillometric arterial blood pressure waveform and points of 

interest therein. 

FIG. 2B illustrates a tonometric arterial blood pressure waveform and points of interest 
therein. 

FIG. 3 is an example display of compliance measurements for analysis according to 
20 the present invention. 

FIG. 4A is a graph of shear force. 

FIG. 4B is a schematic diagram of the forces on an arterial cross section. 
FIG. 4C is a schematic diagram of the forces on an arterial longitudinal section. 
FIG. 4D is a schematic diagram of the forces on an arterial cross section in an arm. 
25 FIG. 5 is an overview of an example process for waveform analysis according to the 
present invention. 

FIG. 6 is a more detailed flowchart of an example process for pulse contour analysis 
according to the present invention. 
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Description of Preferred Embodiments 

In the following detailed description of the preferred embodiments, reference is 
made to the accompanying drawings that form a part hereof, and in which are shown by 
way of illustration specific embodiments in which the invention may be practiced. It is 
5 understood that other embodiments may be utilized and structural changes may be made 
. without departing from the scope of the present invention. 

Figure 1A illustrates one embodiment of vascular-compliance-determining system 100 
according to the present invention. System 100 includes an oscillometric sensor 150 and 
a tonometric sensor 140, both coupled to computer 120. Tonometric sensor 140 provides 

10 tonometric signal 149 to computer 120. In one embodiment, tonometric signal 149 is an 
analog signal that is sampled and analog-to-digital converted by computer 120 at a fixed 
sampling rate (e.g., two hundred samples per second) to provide a series of digital values 
representing the pressure measured at the radial artery of patient 99. Oscillometric sensor 
150 includes cuff 158, pressure sensor 153, and pump 154, and generates oscillometric 

15 signal 159. Pump 154 is controlled by computer 120, and provides both inflation and 
deflation functions. In one embodiment, oscillometric signal 159 is an analog signal that 
is sampled and analog-to-digital converted by computer 120 at a fixed sampling rate 
(typically fifty samples per second) to provide a series of digital values representing the 
gauge pressure of the cuff surrounding the brachial artery of patient 99 as the relatively 

20 DC pressure on the cuff 158 is varied by pump 154. In other embodiments, non-fixed 
sampling rates are used. 

Figure IB illustrates details of one embodiment of oscillometric sensor 150 
according to the present invention. Oscillometric sensor 150 includes cuff 158 
surrounding brachial artery 96 within the patient's arm 95, pressure sensor 153 coupled to 

25 the air line connecting pump 1 54 to cuff 158, and generates oscillometric signal 1 59. 

Artery 96 is shown here in a super-systolic condition, wherein cuff 158 has flattened the 
artery 96, with a tapered artery portion 97 under proximal cuff end 151, and another 
tapered portion under the distal cuff end 152. Pump 154 is controlled by computer 120, 
and provides both inflation (i.e., filling cuff 158 with air) and deflation (i.e., allowing air 

30 to exit cuff 158) functions. In one embodiment, inflation is provided by an air pump, and 
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deflation is provided by a air-bleed valve, both within pump 150 and controlled by 
computer 120. 

Figure 1C illustrates schematically details of one embodiment of cuff 158 over 
an artery 96. 

5 Figure 2 A illustrates an oscillometry arterial blood pressure waveform 210 of 

cuff pressure versus time and points of interest therein. At the left edge of the graph, cuff 
158 has been inflated to super-systolic pressure, however a small amount of oscillation 
exists due to shear stresses in the tapered brachial artery under the proximal end of the 
cuff. 

10 Figure 2B illustrates a tonometric arterial blood pressure waveform and points 

of interest therein. 

Figure 3 is an example display of compliance measurements for analysis 
according to the present invention. 

Figures 4A, 4B, 4C, and 4D show schematically the parameters used in 

15 embodiment 3, described below. Figure 4 A is a graph of shear force. Figure 4B is a 
schematic diagram of the forces on an axial cross section. Figure 4C is a schematic 
diagram of the forces on an axial longitudinal section. Figure 4D is a schematic diagram 
of the forces on an axial cross section in an arm. 

Figure 5 is an overview of an example process for waveform analysis according 

20 to the present invention. Method 500 includes (a) identifying 510 the physiological 
parameter to be quantitatively monitored and estimated, (b) measuring 520 an 
oscillometric signal and a tonometric physiological signal, which signals are 
quantitatively dependent on a particular value for the physiological parameter, (c) 
obtaining 530 a sequence of values that are based on the oscillometric signal and the 

25 tonometric signal, (d) receiving and processing 540 the input signals within the computer 
system to convert the sequence of values to an output signal corresponding to the 
particular value of the physiological parameter, e.g., vascular compliance. 

Figure 6 is a more detailed flowchart of an example process for pulse contour 
30 analysis according to the present invention. This method includes 
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(a) identifying the physiological parameter to be quantitatively monitored and estimated; 

(b) measuring an oscillometric signal and a tonometric physiological signal, which 

signals are quantitatively dependent on a particular value for the physiological 
5 parameter; 

(c) obtaining a sequence of values that are based on the oscillometric signal and the 

tonometric signal; 

10 (d) receiving the sequence of values as input signals to a computer system; and 

(e) processing the input signals within the computer system to convert the sequence of 

values to an output signal corresponding to the particular value of the 
physiological parameter. 

15 

(f) using an oscillometric signal to calibrate tonometric pressure signals in a contralateral 

arterial site. 

In some embodiments, a calibrated radial pressure waveform P r (t) is derived from the 
20 tonometric signal S r (t) as follows: 

P r (t)=(l/ar)(S r (t)-b r )+p 

where a, = ( S r (t D )-S r (t M ) ) / ( DBP-MBP ) , 
25 b r - S r (t M ) - a r MBP , and 

p= gh are calibration factors, and where 
= density of blood, 
g = acceleration to gravity, 

h = height difference between the oscillometric and the tonometric measurement 
30 sites, and is zero if the patient is supine, 

MBP is oscillometric mean arterial blood pressure measured at time t M , and 
DBP is oscillometric diastolic blood pressure measured at time t D . 

(g) calculating a first compliance value based on the calibrated radial pressure waveform; 

35 

(h) estimating end-effects of the oscillometric signal; and 

(i) correcting the first compliance value using the estimated end effects. 

40 As used herein, arterial compliance is defined as the absolute increase in arterial 

lumen volume for a given increase in arterial lumen transmural pressure, or 
Compliance = dV/dP T R. 

The static pressure- volume relationship in the arterial lumen is the result of a 
balance of forces. Forces working to reduce lumen size are atmospheric pressure around 



7 



It . 
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the arm, tensions in the tissues surrounding the artery, and tensions in the walls of the 
artery. Forces working to open the artery are lumen pressure, tensions in the arterial wall 
working to maintain structure and shape, tethering forces between the walls of the artery 
and the tissues surrounding the artery, and shear forces that occur when the expansion or 

5 compression of one part of the artery pulls on another part. In vivo, the pressure wave 
travels at many meters per second, and so an artery has the appearance of increasing and 
decreasing simultaneously along its entire length. This, along with the narrow pressure 
range and relatively flat compliance range that occurs under physiological conditions, 
means that shear forces play a negligible role. However, when a short segment of artery 

10 is compressed underneath a blood pressure cuff, shear stresses occur in the arterial walls 
beneath both ends of the cuff. The effects of these shear forces are for the portions of the 
artery away from the cuff to pull open the portions just beneath the cuff ends. Thus, even 
under complete occlusion by the cuff, when the cuff is far above systolic pressure, 
pressure pulsations from the heart are able to expand the tapered region beneath the 

1 5 proximal end of the cuff because the systolic pressure is augmented by shear stresses. 
The result is that oscillometric signals persist at super-systolic cuff pressures. Using a 
blood pressure cuff to measure arterial compliance thus produces compliance estimates 
that are too high because they include shear stresses that are not present without the cuff 
Without cuff compression, compliance is primarily a function of arterial wall stresses and 

20 anatomical dimensions. 

Cuff oscillometry is currently a widespread technique for obtaining the values 
of one, two, or three blood pressure points during the cardiac cycle. The technique is 
most commonly applied to the brachial artery. In this technique, arterial occlusion is 
obtained by placing a fluid-filled (usually air-filled) compliant chamber such as a cuff or 

25 patch around an artery. The pressure in the fluid is mechanically manipulated and 
measured. Pressure pulsations in the artery create volume displacements which are 
captured by the cuff. When the volume of the cuff changes, the pressure of the cuff fluid 
changes. The pressure is measured with a pressure transducer. As the amount of fluid in 
the cuff is varied between a point producing complete occlusion of the artery, and a point 

30 where the cuff and artery become decoupled, a characteristic oscillometric pressure 
waveform is recorded from the cuff fluid. At high cuff pressures (near arterial 



occlusion), the artery is in a region of low compliance with the rigid cuff limiting the 
distensibility of the artery, and the physiologic pressure pulsations produce low amplitude 
oscillometric signals. As the transmural pressure in the artery approaches the region of 
maximal compliance, high-amplitude oscillometric signals are produced. As transmural 
5 pressures approach their physiologic values, another low compliance region is reached. 
At this point, the low pressures in the cuff make the cuff very compliant. The 
combination of low arterial compliance and high cuff compliance again produces low- 
amplitude oscillometric signals. Pattern-recognition techniques are then used to estimate 
one or more of the following: systolic pressure, diastolic pressure, and mean pressure. 

1 0 Mean arterial pressure is usually estimated as equal to the cuff pressure near where the 
oscillations are a maximum. 

Cuff oscillometry has been used to estimate elastic properties of the brachial 
artery such as compliance (Brinton et al., 1997; Brinton et al., 1996; Chio et al., 1989; 
Chio et al., 1998; Shimazu et al, 1989; see reference list below). The theoretical basis for 

1 5 this capability is that the necessary volume changes are obtained continuously from the 
oscillometric signal while pressure changes are obtained by identifying patterns within 
segments of the signal. Unfortunately, continuous pressures are not obtained. One way 
this may be overcome is to assume a parametric model for the volume pressure 
relationship and then fit this model using the few pressures that can be obtained from the 

20 cuff. Most recent efforts have been to obtain better arterial pressure estimates from the 
cuff pressure, and an empirical compliance curve is assumed. See U.S. Patent No. 
5,579,778 issued Dec. 3, 1996, entitled "Method and apparatus for producing 
thermodilution cardiac output measurements utilizing a neural network" by P. D. Baker et 
al, and U.S. Patent No. 5,339,818 issued Aug. 23, 1994, entitled "Method for 

25 determining blood pressure utilizing a neural network" by P. D. Baker et al. 

The present invention shows how the few known pressure points from the 
oscillometric signal can be coupled with a tonometric signal to estimate a person's 
compliance curve at both physiological and sub-physiological pressures. The latter may 
be used, for instance, to determine the amount of blood-pressure reduction needed in an 

30 individual to reduce their compliance to normal levels. 
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A balance of forces on the artery produces the following relationship {SEE 
Derivations 3}: 

P + P sh = P c + P c ,ex + Pt 

where P is lumen pressure, P S h included all pressures to shear forces, P c is the cuff 
5 pressure exerted on the arterial lumen, P c , ex is the cuff pressure exerted on the tissues 
outside the lumen, and P t includes all circumferential stresses. Transmural pressure is 
defined as 
Ptr = P-Pc 

Combining these relations one gets 

10 P T R = Pc,ex + Pt-Psh 

Under in vivo conditions, or with a very long cuff, P S h = 0. It can be shown that {see 
Derivation 6} P sh = Pi d 2 A/ds 2 + p 2 1/A (dA/ds) 2 

where A is lumen cross-sectional area, s is the position along a line through the center of 
15 the lumen under the cuff, and pi and p 2 are functions of shear moduli and geometry. The 
third embodiment proposes estimating P S h to produce a corrected estimate of transmural 
pressure: 

PtR,c = PtR + Psh 

20 In a First Embodiment of the invention, both oscillometric and tonometric 

signals are measured (in one embodiment, both signals are measured substantially 
simultaneously to obtain a series of value-pairs in order to obtain values that are 
correlated to one another), and a calibration is performed. In one such embodiment, the 
oscillometric signal is measured by a pressure sensor coupled to a pneumatic cuff whose 

25 average (i.e., constant, or DC) pressure is manipulated by the control computer (wherein 
the pressure sensor measures the time-varying (i.e., AC) portion), and the tonometric 
signal is obtained by a microphone sensor and/or sensor holder and wrist brace as 
described in the above mentioned U.S. patent applications Serial No. 09/045,018, Serial 
No. 09/394,301, and Serial No. 09/045,449 (each of which is incorporated herein by 

30 reference). In some embodiments, the tonometric signal is measured at the contralateral 
radial artery site (e.g., the left wrist) from the cuff measurement site (e.g., the upper right 
arm). For the Radial Signal (S r ) as measured, 
S r = a r P r + b r 

where P r is radial pressure, and ^ and b r are calibration constants (multiplicative and 
35 additive, respectively), and the corresponding time-varying waveform is 
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S r (t) = a r P r (t) + b r . 

In some embodiments, two pressure points are needed for calibration. In one 
such embodiment, MBP (mean arterial blood pressure) and DBP (diastolic blood 
pressure) are chosen, and respective corresponding times tM and to (shifted to the 
5 appropriate time within a cardiac cycle, e.g., t D is shifted to the nearest pulse minimum) 
are used: 

S r (t M ) = a r MBP + b r 

S r (t D ) = ar DBP + b r 

a r = ( S r (t D )-S r (t M ) ) / ( DBP-MBP ) 
10 b r =S r (t M )-a r MBP 

The radial artery is approximately six inches below the brachial artery. This creates a 

hydrostatic pressure head that can be accommodated by a further pressure head correction 

factor of p= gh, where 
= density of blood (= 1.03 g/cm 3 ) 
15 g = acceleration to gravity (= 980 cm/sec ) 

h = 6" 15 cm 

thus p =1.03 g/cm 3 980 cm/sec 2 15 cm = 15141 g/cm sec 2 

1 1 mm Hg if sitting (or 0 if supine) 
and the calibrated radial pressure waveform P r (t) is derived from the Radial Signal S r (t) 
20 as follows: 

P r (tMl/aO(Sr(t)-b r )+p 

In another such embodiment, the two pressure points chosen for calibration are 
25 MBP (mean arterial blood pressure) and SBP (systolic blood pressure), and respective 
corresponding times t M and ts are used: 
S r (t M ) = a r MBP+b r 
S r (t s ) = a r SBP + b r 

^ = ( S r (t s )-S r (t M ) ) / ( SBP-MBP ) and 
30 b r = S r (t M ) - ^ MBP. 

In yet another such embodiment, the two pressure points chosen for calibration 
are SBP (systolic blood pressure) and DBP (diastolic blood pressure), and respective 
corresponding times ts and to are used: 
35 S T (t s ) = ar SBP+br 

Sr(tD) = ar DBP + br 
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ar = ( SXts)-S r (tD) ) / ( SBP-DBP ) and 
b r = SAt D ) - a, DBP. 

In a Second Embodiment of the invention, a first compliance value (herein also 
5 called "an uncorrected compliance value") is calculated as follows: 
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u = dP m dP TR >o 



V - V T0TAL - V c - Vjjssue 



dV dV r 



dP TR dP TR 



1 5 In some embodiments, it is assumed that the cuff is filled with an ideal gas. Temperature 
can be eliminated from the ideal gas state equation by assuming adiabatic conditions. 
According to other researchers (for example, Shimazu, Kawarada): 



(Wl= constant = « )5 = 1.4 => V c = f-^_l 



PatmVq _ „ 

— = n o 



20 ^ATM 

(Ideal gas law, before inflation) 

" a - V o KA ATM 

n o 

V 0 = residual volume in pump system and hose 
R = universal gas constant 
25 =8.314 in SI units 

Tatm = Ambient temperature 
P ATM = Ambient pressure 
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The computer system controls a cuff7pump system that inflates and deflates the 
cuff. To keep track of the amount of air in the cuff (lie), one needs to know the inflation 
and deflation rates for the cufFpump system. 



[lt\ fl = fwi(t)=>Dc(t) = n ° + ^ f infi(ndt' 



(inflation) 

dn, 



= fden(t) => n c (t) = nmax + lofdeflO')^ 

dett 



dt 

(deflation) 

n m ax = n 0 + J 0 Tw 'f infl (f)dt 

10 

Examples 

f infl (t) = k^n c (t) = n 0 + kt 
(inflation) 

fdeflW = kn c x> n c (t) = n 0 +(n max -n 0 )e" rt 

15 (deflation) 

Constant r can be obtained for a particular oscillometric device as 

1 0.01n o 
r = In — ,t 01 x 

t.01 n ma.x ~ n o 

time to 99% complete deflation. 



20 In some of its embodiments, the Second Embodiment uses a method employing 

the following procedures: 

1 . Apply cuff oscillometric system to left upper arm 

2. Apply tonometer system to right radial artery. 

3. Simultaneously, collect oscillometric and tonometer signals. 

25 4. Using the oscillometric algorithm, identify MBP and SBP from oscillometric 

signal. These pressure points coincide with time points tM and ts, respectively. 
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The process 50 begins with digitizing an analog tonometric blood pressure 
waveform (52), preferably at 200 Hz, with a 16-bit resolution. According to one example 
embodiment, waveforms are collected for a thirty (30) second duration. These data are 
stored in RAM 22 or in storage device 24. In addition to tonometric blood-pressure- 
5 waveform data, it is also necessary to obtain oscillometric cuff pressure data, input at port 
31, which, in one embodiment, is obtained at the same time as tonometric waveform 
acquisition in order to calibrate the waveform in terms of pressure in mm Hg. 

Next, the data is filtered (54), in some embodiments using an eight (8) pole, 25 
Hz low pass Butterworth filter with an attenuation of 6.0206 dB at the corner frequency, 

10 for the purpose of marking individual beats. The initial filtering at 25 Hz provides 
sufficient smoothing of the blood pressure data to enable consistent marking for the 
upstroke (A), peak (B), and end diastolic point (D) for each beat. These locations are 
identified without the use of derivatives of the pressure data. In order to identify the 
dicrotic notch (C) and the left and right scan points bracketing the notch (L-R), termed 

15 the near-notch region (NNR), derivatives of the pressure waveforms are utilized. Digital 
derivatives inherently amplify noise in the data and thus, to minimize the noise in the 
derivatives of the blood pressure waveforms, the portion of the waveform between 
systole and end diastole for each beat is further filtered in some embodiments using an 
eight (8) pole, 15 Hz low pass Butterworth filter with an attenuation of 6.0206 dB at the 

20 corner frequency. 

After the marks have been identified, the digitized waveform data is converted 
to pressure in mm-Hg for analysis using a two-point calibration. In one embodiment, the 
median systolic and median upstroke values are used together with the systolic and 
diastolic measurements of blood pressure obtained from the oscillometric cuff pressure 

25 measurement unit 35 to calibrate the pressure data. In another embodiment the median 
mean waveform value and the median upstroke waveform value together with the mean 
pressure and diastolic pressure determined from the oscillometric cuff pressure 
measurement unit 35 could be utilized to calibrate the data. Using a preferred tonometric 
method (e.g., as in the f 177 patent), capture uncalibrated pressure signals from the radial 

30 artery at a sampling rate sufficient to identify the peaks and nadirs of individual pulses. 
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Using this signal, identify the peaks and nadirs of the individual pulses using a preferred 
algorithm (e.g., as in the 433 patent). 

1 . Using time points from procedure 4, locate corresponding tonometric signal 
values shift to the nearest peak (for t s ), nadir (for t D ), or compute the mean of the current 

5 pulse (for tu), and calibrate using the First Embodiment, but exclude pressure head 
correction factor p. 

2. Using tonometric and oscillometric pressures, P and P c , compute transmural 
pressure Ptr = P - P c at each time point. 

3. Using P c and lie (see Eqn . VP n above) compute V c . 

10 4. Numerically differentiate the data pairs (-V c , Ptr) to obtain C = 

dV = dV c 
dP TR dP TR 

as a function of Ptr. 

5. Plot C(P TR ) and report C(SBP), C(DBP), C(120), C(80), and pressure at Cmax- 

15 Compliance Range: {C(DBP)} to {C(SBP)} 

Compliance Range at NormaUzed Pressure: {C(80)} to {C(120)} 

Mean Compliance: 

^ SBP 

For example, with SBP=120, DBP=80--> 

j 120 

25 Maximum C: {Cm} 
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Third Embodiment - Correction for Shear 

P = arterial lumen pressure 

P c = cuff pressure 

(7 t = tissue circumferential stress 

O w = wall (of artery) circumferential stress 

a S ht = tissue shear stress 

O shw = wall shear stress 

r = lumen radius 

h w = wall thickness 

h t = tissue thickness 

A x = cross-sectional areas 

The following is Derivations 1-6. 

Please refer to drawings of Figures 4A, 4B, 4C, and 4D. 
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The present invention addresses estimating the systemic vascular resistance 
element of the 3 rd -order modified windkessel model. Three different solutions to this 
problem are presented. 

5 

Some Basic Fluid Dynamic Preliminaries for the Present Invention: 

Two independent state variables are used to determine the properties of an 
arterial system-pressure and flow. The relationships between these variables are 
quantities that are frequently of interest when describing a person's arterial system. For 

1 0 example, the ratio of flow to the time rate of pressure change is arterial compliance, the 
steady-state ratio of mean pressure to mean flow gives vascular resistance, and the 
pressure rate of change and flow rate of change are related through the inertance. 
Because blood flows in an axial direction, it is often useful to express flow as the product 
of cross-sectional area and blood velocity. 

15 Because pressure and flow are both needed to describe system properties, and 

since it is the unknown dependence between these two quantities that is sought by the 
present invention for each individual who is assessed, a measurement must be 
independently made of both quantities. In the 3 rd -order modified windkessel model, as 
described by Cohn & Finkelstein, the problem of obtaining a flow measurement is 

20 simplified by only considering the part of the system where net flow from the heart is 
zero. Since, in addition, systemic vascular resistance is assumed to be constant 
throughout the heart cycle, we need only a measure of mean flow, otherwise known as 
cardiac output. Cohn & Finkelstein obtain an independent measurement of cardiac 
output using an empirical model that is based upon timing information from the pressure 

25 wave, and demographic information about the patient. The present invention provides 
additional approaches that are used, in various embodiments, to obtain an estimate of 
cardiac output. 

In developing a model for cardiac output for use in the 3 rd -order modified 
windkessel, there are three goals. First, the present invention is concerned only with 
30 prediction and not interpretation. That is, the model may be treated as a black box with 
any variables found to be useful treated as inputs and with cardiac output as the black box 
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output. Second, the model needs to improve upon a model that attributes the same 
cardiac output (e.g. 5 liters per minute) to all individuals. Even with such a crude model 
of cardiac output, we know that the windkessel provides discriminating information about 
arterial compliance. Therefore, improved models of cardiac output may be considered as 
5 methods of fine-tuning the modified windkessel. Finally, the present invention makes 
use of more of the available information that might be useful. This includes the blood 
pressure and patient demographic information. The conventional Cohn/Finkelstein 
cardiac output model (see, e.g., US Patent No. 5,241,966 and Cohn J.N., Finkelstein S., 
McVeigh G., Morgan D., LeMay L., Robinson J., Mock J., "Noninvasive pulse wave 

10 analysis for the early detection of vascular disease," Hypertension. 26(3):503-8, 1995 
Sep.) meets these three criteria, and next various embodiments are described using three 
additional approaches. 

Before detailing the models, it is useful to recall a deficiency in the 3 rd -order 4- 
element modified windkessel model. For mathematical clarity, the lumped electrical 

15 analog is usually illustrated, rather than the mechanical analog of the model. This is 
shown in Figure 1 . 

A mass balance through the system leads to the following equations: 
Equation (1): io = ici+iL 
Equation (2): iL = ic2+iR 

20 Equation (3): i C i=CidPi/dt 
Equation (4): icr=C 2 dP 2 /dt 
Equation (5): L dii/dt = P 2 -Pi 
Equation (6): P 2 =R iR 

25 where the i's are flows through respective elements. After some minor algebraic 
manipulations we obtain the following 3 rd -order ordinary differential equation: 

Equation (7): Ri 0 = { 1 + (RCi+RC 2 )D - (RC,L/R)D 2 - (RC, L/R RC 2 )D 3 }P 2 

30 where D is the differential operator with respect to time. Setting i 0 =0 (for diastole) and 
using Laplace transforms we easily obtain an expression for P 2 as a function directly of 
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time, and indirectly of RQ, RC 2 , and L/R. We then use curve-fitting to estimate RCi, 
RC 2 , and L/R from a sample of measured P 2 data. Herein lies the problem--we are unable 
to estimate Ci, C 2 , or L without a separate estimate of R. 

To obtain an estimate of R, we again use equation (6). Throughout we have assumed that 
R is constant, so integrating both sides we obtain JP 2 dt = R J i R dt. We have P 2 data, but 
we don't know iR, so we rearrange equations (l)-(2) to obtain i R = (i 0 - ici) - ic2- 
Integrating, we get \ i R dt = j i 0 dt - ( \ i C idt + j i C2 dt). Substituting equations (3)-(4) we 
get Ji C idt+ ii C2 dt= jCidP,/dtdt+ jC 2 dP 2 /dt dt = C\ fdPj + C 2 JdP 2 = 
Ci(Pi - Pi jt =to) + C 2 (P 2 - P 2 ,t=to)- Since Pi, like P 2 , is a periodic function at steady state, it 
reaches nearly the same pressure value at any two points separated by a pulse duration. 
We get C,(P, ,t=to+T - Pi,t=to) + C 2 (P 2>t =to+T - P2,t=to) - 0, so that finally, 

(8) Equation (8) R- JP 2 dt/ \ i 0 dt, 

where the integrals are over a full cardiac cycle. The reason this expression for R is 
likely to be useful is that 

(9) Equation (9) \ i 0 dt = CO/HR = stroke volume. 

Heart rate (HR) we can easily obtain from the pressure wave, so we can choose to 
estimate either stroke volume, or cardiac output. 

The Cohn/Finkelstein model for stroke volume is a direct linear parametric model. Here 
we offer three additional categories of models — an indirect linear parametric model, an 
indirect nonlinear nonparametric model, and a direct nonlinear nonparametric model. 

Indirect Linear Parametric Model 

We call this approach linear parametric because it uses ARMAX (AutoRegressive 
Moving Average with exogenous input) models [1]. It is indirect because instead of 
estimating stroke volume directly, ascending aortic cross-sectional area and blood 
velocity are modeled and combined to estimate stroke volume. 
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' Two ARMAX models are used, one to compute a transfer function from radial arterial 
pressure to ascending aortic velocity, another to estimate ascending aortic cross-sectional 
area from radial arterial pressure. A brief review of ARMAX modeling is next presented. 

5 An ARMAX model for y(t) as a function of x(t) is given by 

Equation (10): y[n] = E p =i, P a[p]y[n-p] + E q = 0> Qb[q]x[n-q] + w[n] + E F i, R c[r]w[n-r] 

where w is the unknown error, and the vectors a, b, and c are the model parameters. This 
1 0 equation can be rewritten in matrix form as 

Equation (11): Y = Fa + WB 

where the matrix F consists of columns of shifted vectors of y and x, a contains the V 
15 and 'b' parameters, and 6 contains the 'c' parameters. A common way of obtaining 
estimates of the a, b, and c parameters is to use an iterative two-step ordinary least 
squares (OLS) algorithm followed by a generalized least squares (GLS) algorithm until 
convergence is reached [2] . 

20 Taking the Z-transform of the above expression for y[n] and rearranging, we obtain 

Equation (12): y(z) = (B/A) x(z) + (C/A)w(z) 

where B is a polynomial in z" 1 with coefficients 'b\ A is a polynomial in z" 1 with 
25 coefficients 'a', and C is a polynomial in z" 1 with coefficients *c\ Taking the inverse Z- 
transform gives 

Equation (13):y[n] = Ej= UN h v [n-i]x[i] + E j=i jN h w [n-i]w[i]. 
30 Estimating y[n] by its expected values (i.e. neglecting error) gives the desired result 
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Equation (14): y[n] = E \ h v [n-i]x[i]. 

ARMAX modeling is straightforward except for the selection of the proper order 
(P, Q, R). A common way of comparing models of different orders is with a bias- 
5 adjusted maximum likelihood estimate such as AIC. Several candidate models are fit and 
the one with the lowest AIC is chosen. 

Allen & Murray used ARMAX models to predict photoelectric 
phlethysmographic blood volumes from isovolumic calibrated photoelectric finger 
pressures [3]. They used a first-order two-parameter autoregressive (ARX) model (P, Q, 
10 R) = (1, 1,0) with b[0]=0. The ability of these low-order models to predict volumes from 
pressures was modest, but suggests that reasonable models might be obtainable with only 
a few parameters. Allen & Murray related finger volumes to finger pressures, while our 
needs are to relate aortic variables to radial arterial pressure. Thus the transfer functions 
that we compute implicitly combine a local pressure to velocity (or area) transfer function 
1 5 with an inter-positional transfer function (aorta to radial artery). 
We obtain stroke volume as 

Equation (15): SV = |flow(t)dt = jA(t)v(t)dt E t n= i, N A[n]v[n] 

20 where A(t) is ascending aortic cross-sectional area as a function of time, v(t) is ascending 
aortic blood velocity as a function of time, A[n] and v[n] are the corresponding sampled 
data values, t is the sampling interval, and N is the number of data points in a single 
cardiac cycle. Here the integral is approximated with a simple sum, but any appropriate 
numerical integration could be used to obtain higher precision (e.g. a high-order Newton- 

25 Cotes). A[n] and v[n] are obtained as 

Equation (16): A[n] = E j h A [n-i]P 2 [i] 
Equation (1 7): v[n] = E i h v [n-i]P 2 [i] 
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where P 2 [n] is the sampled pressure at the radial artery. The transfer functions h A and h v 
are obtained using two separate ARMAX models as described above. The model search 
may be restricted to the family (P, Q, 1). 

5 The following procedure is therefore proposed for any particular patient group of 

interest: 

1 . Using a Doppler ultrasound (US) probe, obtain real-time ascending aortic cross- 

sectional areas and blood velocities while simultaneously collecting radial arterial 
10 tonometric pressures. Obtain an oscillometric cuff mean and diastolic blood 

pressure reading immediately before the collection. Obtain several cardiac cycles 
of data for each of several human subjects. 

2. For each subject, estimate the transfer function parameters (a's, b's, and c's) using the 

two-step ARMAX algorithm described above and selecting the model with the 
1 5 lowest AIC. Obtain both the pressure to velocity transfer function and the 

pressure to area transfer functions. 

3. Average the transfer function parameters for the several beats and for all individuals. 

4. Using the averaged transfer function parameters, compute the two transfer functions 

(h v [n], and h A [n]). 

20 5. In all future windkessel analyses, estimate SV using Equations (15), (16), and (17). 

Indirect Nonlinear Nonparametric Model 

This method is called indirect because, again, area and velocity are estimated 
from pressure. It is called nonlinear nonparametric because a neural network is used. 

25 This method is similar to the indirect linear parametric model, with two main differences. 
First, we train a neural network rather than an ARMAX model. Second, because a neural 
network is being used, additional parameters are easily incorporated. Similar to the 
ARMAX model, we need to decide upon a model order (Q). In addition, we need to 
decide which additional variables to include. As with the present PulseWave® 

30 implementation, we use gender, age, height, and weight. Data are collected with Doppler 
ultrasound and arterial tonometry as described above. Any network activation type can 
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be used, but for demonstration purposes, we will describe a logistic two-layer multiple- 
layered perceptron (MLP) model. The number of inputs is equal to 4 + Q + 1 . There are 
two outputs (velocity and area). The number of nodes between the two middle layers 
must be high enough to describe the input-output map, but not so high as to make training 
5 or prediction computationally impractical. We always use the same number of nodes in 
both layers. The number of nodes is between one and two times the number of inputs. 
Training proceeds with cross-validation until leave-one-out cross validation sum of 
squared residuals reaches a minimum (rather than until actual sum of squared residuals 
reaches a minimum). It is because we are unconcerned with the parameter values of the 

10 neural network that this approach is called "nonparametric". 

Any nonlinear parameter estimation algorithm may be used, e.g. conjugate 
gradients, backpropagation, Newton or quasi-Newton methods. Training stops at 
minimum leave-one-out cross-validation sum-of-squared residuals. The cross-validation 
sum-of-squares is also used to compare models of different orders. For each beat, we 

15 train the two outputs (v[n] and A[n]) to the corresponding pressure, P 2 [n], as well as the 
previous Q P 2 values, P 2 [n-i] for i = 1 . . .Q. Age, height, weight, and sex are also input. 
The model can be written as 

Eq. (18): (v[n], A[n]) = NN(age, height, weight, sex, P 2 [n], P 2 [n-1], . . ., P 2 [n-Q]) 

20 

Therefore the following procedure is used, in one embodiment, for a particular group of 
interest: 

1 . Collect data as in the indirect linear parametric model above. Because of the 

25 generality of neural network models, it may be necessary to obtain data from 50 

to 100 individuals. 

2. Using the method described above, train and select a neural network configuration 

that predicts v[n] and A[n] from P 2 [n], age, height, weight, and sex (Equation 
(18)). 

30 3. In all future windkessel analyses, estimate SV using Equations (15) and (18). 
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Direct Nonlinear Nonparametric Model 

This method is called direct because S V is computed directly. It is called 
nonlinear nonparametric because a neural network model is employed. 
This approach is similar to the indirect nonlinear nonparametric model except that the 
5 neural network output is SV, and a fixed set of pressure values are used as inputs. The 
model is expressed in the following way: 

Equation (19): SV = NN(age, height, weight, sex, P 2 [0], P 2 [l], • . ., P 2 [tf]) 

Unlike with the indirect nonlinear approach, this approach is dependent upon the 
start and end times used for the pressure inputs. We use the foot of each beat for the first 
10 pressure point. Various end-points can be tried to obtain the lowest prediction error. The 
spacing between samples can also be adjusted to obtain simpler models with fewer 
inputs. Training should otherwise proceed as described above for the indirect nonlinear 
nonparametric model. 

The following procedure is used, in one embodiment, for a particular group of 
15 interest: 

1. Using a Fick dilution method, estimate several cardiac outputs for an individual. 

Simultaneously, collect tonometric pressure recordings at the radial artery that are 
calibrated to brachial oscillometric cuff blood pressure estimates. 
20 2. Using the method described above, train and select a neural network configuration 
that predicts cardiac output from P 2 [n], age, height, weight, and sex (Equation 
(19)). 

3. In future windkessel analyses, estimate cardiac output using Equation (19). 
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Conclusion 



One aspect of the present invention provides a method for computerized 
calculation of a variable physiological parameter of a patient. The first method includes 
(a) identifying the physiological parameter to be quantitatively monitored and estimated, 
30 (b) measuring an oscillometric signal and a tonometric physiological signal, which 

signals are quantitatively dependent on a particular value for the physiological parameter, 



24 



1 » 

(c) obtaining a sequence of values that are based on the oscillometric signal and the 
tonometric signal, (d) receiving the sequence of values as input signals to a computer 
system, and (e) processing the input signals within the computer system to convert the 
sequence of values to an output signal corresponding to the particular value of the 

5 physiological parameter. 

In some embodiments, the first method further includes (f) using an oscillometric 
signal to calibrate tonometric pressure signals in a contralateral arterial site. 
In some embodiments, a calibrated radial pressure waveform P r (t) is derived from the 
tonometric signal S r (t) as follows: 

10 P r (t)=(l/ar)(S r (t)-b r )+p) 

where a, = ( S r (t D )-S r (t M ) ) / ( DBP-MBP ) , 
b r =S r (t M )-a T MBP,and 
p= gh are calibration factors, and where 
= density of blood, 

15 g = acceleration to gravity, 

h = height difference between the oscillometric and the tonometric measurement sites, 

and is zero if the patient is supine, 
MBP is oscillometric mean arterial blood pressure measured at time t M , and DBP is 
oscillometric diastolic blood pressure measured at time t D . 

20 In some embodiments, the first method further includes (g) calculating a first 

compliance value based on the calibrated radial pressure waveform, (h) estimating end- 
effects of the oscillometric signal, and (i) correcting the first compliance value using the 
estimated end effects. 

In some embodiments, the processing function (e) includes estimating a first 

25 compliance value using a compliance pressure curve. 

Some embodiments of the first method further include (j) using a tonometric signal to 
calibrate oscillometric pressure signals in a contralateral arterial site. 
Some embodiments of the first method further include (h) estimating end-effects of 
oscillometric sensor apparatus on the oscillometric signal. 

30 Another aspect of the present invention provides a method for computerized 

calculation of a variable physiological parameter of a patient. This second method 
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includes(a) identifying the physiological parameter to be quantitatively monitored and 
estimated, (b) coupling at least a first sensor and a second sensor to the patient, the first 
sensor and the second sensor each being responsive to monitor different time- varying 
physiological waveforms, which waveforms are quantitatively dependent on a particular 
5 value for the physiological parameter, (c) obtaining a time-correlated dual sequence of 
digital values that are based on the waveforms monitored by the first and second sensors, 
(d) receiving the sequence of digital values as input signals to a computer system, and (e) 
processing the input signals within the computer system to convert the time-correlated 
dual sequence of digital values to an output signal corresponding to a value of the 
1 0 physiological parameter. 

In some embodiments of the second method, the physiological parameter is 
vascular compliance, the first sensor monitors an oscillometric waveform, and the second 
sensor monitors a tonometric waveform. 

In some embodiments of the second method, the physiological parameter is vascular 

15 compliance, the first sensor monitors an oscillometric waveform derived from 
oscillometric signals of a brachial artery site, and the second sensor monitors a 
tonometric waveform derived from tonometric signals of a contralateral radial artery site. 

Yet another aspect of the present invention provides a system for computerized 
calculation of a variable physiological parameter of a patient. The system includes a first 

20 and a second sensor. The first sensor measures an oscillometric physiological signal that 
is quantitatively dependent on a particular value for the physiological parameter. The 
second sensor measures a tonometric physiological signal that is quantitatively dependent 
on the particular value for the physiological parameter. A first analog-to-digital 
converter, operatively coupled to the first sensor, generates a first sequence of digital 

25 values that are based on the oscillometric signal. A second analog-to-digital converter, 
operatively coupled to the second sensor, generates a second sequence of digital values 
that are based on the tonometric signal. A computer system, operatively coupled to the 
first and second analog-to-digital converters, processes the first and second sequences of 
values to generate an output signal corresponding to the particular value of the 

30 physiological parameter. 
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In some embodiments, the first sensor senses the oscillometric signal from one 
side of a patient, the second sensor senses the tonometric signal from a contralateral 
arterial site, and the computer uses the oscillometric signal to calibrate tonometric 
pressure signals in the contralateral arterial site. 
5 In some embodiments, the computer derives a calibrated radial pressure 

waveform P r (t) from the tonometric signal S r (t) as follows: 
P r (t)=(l/ar)(S r (t)-b r )+p 
where a r = ( Sr(t D )-S r (tM) ) / ( DBP-MBP ) , 
b r =Sr(tM)-a r MBP,and 
10 p= gh are calibration factors, and where 
= density of blood, 
g = acceleration to gravity, 

h = height difference between the oscillometric and the tonometric 
measurement sites, and is zero if the patient is supine, 
1 5 MBP is oscillometric mean arterial blood pressure measured at time t M , & 

DBP is oscillometric diastolic blood pressure measured at time t D . 
In some embodiments, the computer system further calculates a first compliance value 
based on the calibrated radial pressure waveform, estimates end-effects of the 
oscillometric signal, and corrects the first compliance value using the estimated end 
20 effects. 

In some embodiments, the computer system further estimates a first compliance 
value using a compliance pressure curve. 

In some embodiments, the computer system further uses a tonometric signal to 
calibrate oscillometric pressure signals in a contralateral arterial site. 
25 In some embodiments, the computer system further estimates end-effects of oscillometric 
sensor apparatus on the oscillometric signal. 
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Another aspect of the invention is to improve the accuracy of the estimate of 
stroke volume using one of three methods. A first method includes measuring a relevant 
population of persons to determine empirically the h A and h v parameters for Equations 
16 and 17 respectively, and then using Equations 16 and 17 to derive arterial cross 
section area A[n] and velocity v[n] respectively. These A[n] and v[n] parameters are 
then input into Equation 15 to determine an improved estimate of stroke volume. A 
second method further obtains age, height, weight and sex information for each person in 
the population used to obtain data to empirically obtain weighting factors, and Equation 
18 is then used to obtain the A[n] and v[n] parameters, which are then input into Equation 
15 to determine an improved estimate of stroke volume. A third method also obtains age, 
height, weight and sex information for each person in the population used to obtain data 
to empirically obtain weighting factors, and Equation 19 is then used to directly obtain an 
improved estimate of stroke volume. These improved estimates of stroke volume are 
then used in the methods described above to obtain improved estimates of compliance 
values. 

On the following pages are explanations of the derivations 1-6 described above. 

The following is Derivation 1 : 
Please refer to drawings of Figures 4A and 4B. 
. Here we derive a general expression for a shear force. 



(1989). 





Hie following is Derivation 2: 

Using Derivation 1, we derive expressions for the rate of change of shear force along the length of the 
artery. We do so for both the wall and tissue layers. We also define the geometric coefficients, 
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4C and 4D, we wri«e a force batoee, and derive re.ariona.np between 
various pressures in title aim and cuff. 

Force balance: 

F + F shw + F slu =F c +F„+F, 

where 
F = Pr„5t 

F shx =Fjt+#)-F shx {t) 

6tF x = cr x h x & 

and x is either t or w. 
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The following is Derivation 4: 

Here we compute simplified expressions for the terms that appear in the geometric coefficients in 
Derivation 2. 

A 

**w 

dh r k J 

Li. 



dh 



The following is Derivation 5: 

Here we uses Derivations 2 and 4 to estimate the value of the geometric coefficients based on the 
assumptions that the ratio hjr w is always approximately 0.06, and the ratio b t /r w is always approximately 

y * 0.06 =>a w * 0.12, 5l * io => 4 * 0.12, 
« 1 .06,, . V w2 * 0.124, ^ * 20, » 222 
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When we measure C = -7=: — , the volume change is due not only to Ptr , 

dPxR 

but to P s ii as well. We want to know what fraction of the total volume change 

is due only to P^R • Normally we are computing the total volume change in the 

cuff or artery (dV = dV to t). We will find a correction factor to convert to the 

dVjR 
desired -pr — . 
dPjR 

dVTR = dVtot / dV ra /dt i = dV^ 
dP TR dPTR I dVtrt/dt I dPxR 



jA(x,t)dx> V TR = LA(0,t) 



dVtot = I A(x, t) dx 



where L is the cuff length. 
Letting r w refer to ^(O, t), 

dt " f dt 
k ~ dVxR/ dt 
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= 1 t A(0 - l) + £ £ *°» '> 

We can find a model for A(x) by solving the nonlinear ODE from Derivation 3. 
d 2 A 



Psh = 



1 4* ^ A I dx / 1 4jt ' 



dx 2 



or 



J 2 ^ 1 /dA\ 2 
where 

1 /'sh2 fva + faiVt/Vw) 



^ = 777" Psh = 7T~ p 4 "V"^ + p t ~ tr) 



Assume that wall and tissue properties are similar and use Derivation 5. Also 
use /u w - 137 kPa = 1027.6 mmHg (from Deng et. al., American Journal 
of Physiology. 266(1 Pt 2):H1-10, 1994 Jan). We get 

/' 1 

Ht « AAf ~ 1027.6 mmHg, <r t ■» cr w « P ~ « P> 311(1 

ori « -0.905 « -1, 

0.0566 i/ / 1 \ \ 

" 2 w T02^ «°Ml0.06 + 0^6 P ) " 

« 10~ 3 mmHg -1 (9.18 /> + 0.609 P c ) 
Substituting ari , 

rf 2 ,4 1 /dAy* . tf 2 . . A . 
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Solving for A(x) and centering around Xq = 0, and taking r w to mean r(0, t), 
get 



A(x) = A w cosh 2 (— 77), 



A X - Xq A / 2 
Next we evaluate each of the terms in κ: 

£ A( - x > trn M( ! + cosh (£ *)) - i £ * x )) 

-r- A{0, t) = 2xr w 

A(x, t) dx = y r w {2 Cr H -^) + 3A-A) 



* ^,0dx = i(^) 2 (A-A) 



da 2 2 v p 

where 



.VL=°'.^ = i/2 

A = x " coi 4* £ ) - *■ C0 H X >- £) = I C0S K4 £) 



Substituting back into ; we get 

K 



3H 



! fr w {2(x H -x L ) + 3D 2 -D 1 } + % \ ( r ff (A " ^2) 



Z 2jrr, 



= _L (2 ,,.., + A) + (A .z> l)(l+ (J-^)(^f)} 
= _l_ {2fe ., I + A) + (z ,_ A)(l+ |gg) } 

= 4^{2to-^ + Z) 2 ) + (i> 2 -A)(l + 



A dlog(af 2 ) _ dlog(P sh ) 
= dlog(r w ) dlog(/ w ) 
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~~ 10" 6 mmHg" 1 (387 P c + 5842 />) 

To estimate 7' w we could use 0.25 cm, or choose sigmoid model (like 
arctangent) with pressure as is sometimes used. Here, for example, we use a 

logistic curve ^logistic(.%) = - - j . 

~ (0.5 cm) logisticf — ±^hr — * 

^ \ 20mmHg / 

In one embodiment, the Third Embodiment uses the following procedure: 

1. Calculate P, P c , Ptr> 311(1 c f° r second embodiment. 

2. Calculate r w using Pjr and the sigmoid model. 

3. Calculate 0 using P c and P. 

4. Calculate *S by numerically differentiating logS 2 with respect to logr^ and 
multiplying the result by ^ . 



3JT 



5. Calculate D\ and Z>2 from — . 

6. Calculate f£ from A,Z> 2 > 311(1 s * 

7. Calculate C C01TC ctcd = AT C. 



It is to be understood that the above description is intended to be illustrative, and 
not restrictive. Many other embodiments will be apparent to those of skill in the art upon 
reviewing the above description. The scope of the invention should, therefore, be 
determined with reference to the appended claims, along with the full scope of 
equivalents to which such claims are entitled. 
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